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Abstract. Using the dynamical system approach, properties of cosmological models 
based on the Hofava-Lifshitz gravity are systematically studied. In particular, the 
cosmological phase space of the Hofava-Lifshitz model is characterized. The analysis 
allows to compare some key physical consequences of the imposition (or not) of detailed 
balance. A result of the investigation is that in the detailed balance case one of the 
attractors in the theory corresponds to an oscillatory behavior. Such oscillations can be 
associated to a bouncing universe, as previously described by Brandenberger, and will 
prevent a possible evolution towards a de Sitter universe. Other results obtained show 
that the cosmological models generated by Hofava-Lifshitz gravity without the detailed 
balance assumption have indeed the potential to describe the transition between the 
Friedmann and the dark energy eras. The whole analysis leads to the plausible 
conclusion that a cosmology compatible with the present observations of the universe 
can be achieved only if the detailed balance condition is broken. 
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1. Introduction 

Recently, Hofava made a proposal for an ultraviolet completion of general relativity 
(GR) normally referred to as the Hofava-Lifshitz (HL) theory pp. (The name is 
due to Hofava's initial inspiration on the Lifshitz theory in solid state physics.) 
The salient characteristic of the HL proposal is that it seems to be renormalizable, 
at least at the level of power counting. This ultraviolet behavior is obtained by 
introducing irrelevant operators that explicitly break Lorentz invariance but ameliorate 
the ultraviolet divergences. On the other hand, Lorentz invariance is expected to be 
recovered at low energies, as an accidental symmetry of the theory. 

Originally, the HL proposal came with the possibility of imposing or not the so- 
called projectability condition and the detailed balance condition. The first condition is 
related to the space-time dependence of the lapse function, N, which characterizes a 
canonical 3 + 1 decomposition of the metric field g, while the second is a restriction 
on the form of the potential terms which may appear in the Lagrangian that leads to 
simplifications since it reduces the final number of couplings. Notice therefore that we 
have, in principle, four different incarnations of this proposal. 

Since its publication, the HL theory has been the object of an exhaustive research 
regarding its different properties and implications to space-time physics. In particular, 
a lot of attention has been paid to its internal consistency, to how to define the infrared 
limit, its compatibility with GR, and the potential application of the results obtained 
to cosmology. 

Presently — as this article is being written — the consistency status of the theory 
is not completely clear, nor its low energy limit and, hence, how GR is recovered at 
the different regimes. In fact, for the non-projectable version, although the low energy 
limit has been found [2], there seem to remain important problems, regarding strongly 
coupled features that may preclude any type of perturbation approach P, Hj, while 
the corresponding symplectic structure and systems of constraints do not seem to be 
well defined [5J [T], [6] . Also, imposing detailed balance leads to a cosmological constant 
with the wrong sign, from what we should expect to be at odds with cosmological 
observations [HE]. In comparison, the projectable version seems to be less problematic, 
since the above listed problems can in principle be evaded by the non-local form 
of the Hamiltonian constraint [7J. Also, if detail balance is not imposed a richer 
phenomenology seems to appear, where cosmological applications may lead to new 
results in inflation, bouncing cosmology, dark matter, and dark energy (see, for example, 

pisiiiDiiniiiaiiainiiisiiisiirziiiB]). 

At this point it is important to investigate the key aspects of the theory, which may 
help in clarifying the status of the different HL proposals as plausible candidates of a 
quantum theory of gravity. By now, it is clear that the study of the consistency of the 
theory is a complicated problem that will still take a lot of effort and time to ascertain 
and that will involve detailed research of the renormalization properties of the theory 
and, most probably, of its non-perturbative nature (see e.g. [T9l [20]). On the other 
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hand, another important requirement the theory should fulfill is its implication with 
phenomenology and, in particular, with cosmology. In this last regard there are already 
papers on some of the consequences of these new scenarios^; however, no systematic 
study of the possible form of the cosmological relevant solutions of the theory, either for 
the detailed balance case or for a general potential, has appeared yet. 

In the present work, motivated by the above facts, we investigate the space of all 
possible classical solutions of the HL theory that are relevant for cosmology, i.e. we 
study the cosmological phase space of the HL model. To address this involved issue, 
we will borrow techniques from the field of dynamical systems that are frequently used 
on more canonical studies applied to diverse types of cosmologies. As a matter of fact, 
dynamical system techniques in cosmology have been used for many years. In our 
work, we will consider the dynamical system approach to cosmology as formulated by 
Collins and applied by Wainwright and Ellis [22]. Such approach has the advantage 
of relying on dynamical system variables which are directly related to cosmological 
observables — like the matter density parameter — and of being relatively easy to apply 
to very complicated cosmological models. During the past years this method has been 
able to unfold some very interesting properties of Bianchi universes [22], as well as of 
scalar-tensor and higher-order gravity cosmological models (for some detailed examples, 
see [251 1221 122 [26]). 

As we will see, when applied to HL cosmologies these techniques will be able 
to reveal some interesting results. In particular, when detailed balance is imposed, 
although one is able to prove that the corresponding cosmology could — given suitable 
initial conditions — behave like a Friedmann one, it turns out that it will eventually 
evolve towards an oscillatory behavior or either recollapse, leaving no space for a dark- 
energy era. When detailed balance is broken, however, the situation changes: given 
the proper initial conditions and values of the parameters the cosmology undergoes an 
almost Friedmann era(s) to evolve towards a de Sitter era which can be associated to 
dark energy. In this sense the plausible consequence is that a cosmology compatible with 
the present observations of the universe can be achieved only if the detailed balance 
condition is broken. In addition to that, we will show that the information on the 
cosmological relevance of the values and signs of the parameters in the theory helps 
to uncover properties of the general theory. For example, from the relative sign of the 
Hubble and matter terms in the cosmological equations, we will be able to deduce that 
the spin zero modes of the theory can be excluded. 

The paper is organized as follows. In Sect. [2] we present the Lagrangian of the HL 
theory for the two cases of either imposing or not the detailed balance constraint, where 
the minimal potential defined in [9] for the non-detailed-balance case will be used. The 
projectability condition does not affect our analysis, since for the cosmological ansatz 
the shift iV is only time dependent. In Sects. [3] and H] we present the specific analysis for 
the theory with detailed balance and for the one without detailed balance, respectively. 

| See [Ml [21] for some general results of this type based on different assumptions and technics. 
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We will be able to characterize the classical phase space, discussing its structure in 
depth, in particular all the fixed points and their nature as repellers or attractors, in 
each of the cases. Finally, in Sect. [5] we summarize the results obtained, giving some 
perspectives for further work. 



W jl = {9jm,l + 9ml,j ~ 9jl,m) • (3) 



2. Hofava-Lifshitz gravity 

The following is a short overview on the minimal set of definitions needed to understand 
the dynamical system analysis of the forthcoming sections. In the HL theory, the 
dynamical variables are defined to be the laps N, the shift iVj and the space metric 
gij, Latin indices running from to 3. The space-time metric is defined using the ADM 
construction, as 

ds 2 = -N 2 dt 2 + g^dx 1 + JSTdt) (dx j + N j dt) , (1) 

where N % = g^Nj as usual. The action S is written in terms of geometric objects, 
characteristics of the ADM slicing of space-time, like the 3d-covariant derivative Vj, 
the spatial curvature tensor Rijki, and the extrinsic curvature K^. They are defined as 
follows, 

R l jk i = W { jlik - W l m + W m jl W i km - W n lk W\ m , (2) 

where W l ji are the Christoffel symbols (symmetric in the lower indices), given by 
1 
2 

The Ricci tensor is obtained by contracting the first and the third indices 

Rij = g kl R ik ji and R = g ij R^ , (4) 

while the extrinsic curvature is defined as 

Kij = ^{-9H + V *^' + V.-iVi) , (5) 

where the dot stands for time derivative. 

In terms of the above tensor fields, the HL action can be written as 

being the kinetic term universally given by 

Cunetic = a(KijKV - \K 2 ) , (7) 

with a and A playing the role of coupling constants. The potential term is, in principle, 
a generic function of Rijki and V». Here we will work with two different types of 
potentials: first, the detailed balance potential and, second, the potential defined in 
[9] (the SVW case), that corresponds to the more general potential with dimensionless 
kinetic couplings (see the original article for clarifications.). Finally, the matter term 
corresponds to the coupling of the matter fields to gravity. 
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The potential for the detailed balance case is 

Cpotentiai-dztaa = PC^C* + ^ k R u V,R l k + (R i:j R ij + V R 2 + 5R + a, (8) 

where JgC^ = e Vfc(i2j — 1/ARgf) and the Greek letters define the different coupling 
constants. We have closely followed the notation of [8] for simplicity. We stress that 
only three out of the whole set of coupling constants above are independent, due to the 
detailed balance constraint. Then, it is not difficult to see that, if we want the theory to 
make contact with GR in the infrared regime, we have to do the following identifications: 

c , , 1 -2A 

a = 1fl n , A = l, 5 = , <r= — , (9) 

with all other couplings being equal to zero. To get the above results, we also 
have defined a new time coordinate x° = ct, where (c,Gn,A) are the velocity of 
light, Newton's constant, and the cosmological constant, respectively. From this point 
onwards, we will keep c = 1, so that our time coordinate is actually t. 
The potential of the SVW case is 

^potential — SVW 

= g&iRjkSrR* + g 7 RV 2 R + g e R l 3 RlR k + g 5 R(R jk R 3k ) + g^R 3 (10) 
+ g 3 R jk R jk + g 2 R 2 + giR + g , 

where we have modified the original notation to better accommodate our formulation of 
the universal kinetic term (JTj), therefore our coupling constants are not all dimensionless. 
As in the detailed balance case, in order to recover GR at low energies we get similar 
equations to those in Just notice we need to identify 8 g± and o <^ g , while all 
the other couplings in the potential should go to zero. 

At this point, other phenomenological and theoretical considerations may help to 
constraint the range of values that the different couplings should take. For example, 
in [29] it was found that ghost instabilities are present if A G (1/3, 1), the cosmological 
constant is negative for the detailed balance potential, a > 0, etc. Here, we will 
constraint as little as possible the different ranges of values on each coupling constant 
to see how much information comes out of the dynamical system approach. Then, we 
will add this information to the constraints arising from other considerations, to finally 
obtain the most promising form of the potential. In particular, we will take A different 
from 1/3 (corresponding to the scale invariant case) as the only limitation on its range. 

To start our studies on cosmological applications, we have to consider a FLRW 
ansatz on the 4D metric. Due to the homogeneity and simple time dependence of this 
type of ansatz, only a subset of the coupling constants plays a role in the dynamics, in 
any case. Also, as we already said, is a function of time only and, therefore, all issues 
related to the projectability condition turn out to be irrelevant^. Then, we set, 

N^N(t), Ni^O, g .. —mty-yy , (11) 

§ Recall that we are interested in the classical phase space and its structure. On the other hand, the 
dynamics of fields propagating in the above solutions are determined by the specific incarnation of the 
HL theory chosen, where projectability is important. 
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where 7^ is a maximally symmetric 3D metric, of constant curvature R = 6k, with 
k = (-1,0,1). 

The inclusion of matter content in the theory in its general form has not been worked 
out yet. There are various discussions and suggestions (see for example [8]), but this 
is still an open subject. An interesting approach is based on adding some new degrees 
of freedom to the framework, so that the theory gains in diffeomorphism invariance. 
Then, minimal coupling to matter is invoked, as in GR, although new couplings to 
the extra degrees of freedom could also be introduced (see [U [3QI [31] for this type of 
constructions). Here, we will follow previous discussions where a very general approach 
is taken. We will add to the gravity field equations a cosmological stress-energy tensor, 
such that in the low-energy limit we recover the usual GR formulation. This tensor is a 
hydrodynamical approximation with two quantities, p,p, corresponding to the density 
and pressure, and where both matter fields are related by the usual equation of state 
p = wp. Since one of our goals is to investigate the relation between HL and dark 
energy (cosmic acceleration), we will only consider here w > 0, so that no dark energy 
is introduced by hand in the model. 



3. Hofava-Lifshitz cosmology: detailed balance case 

From the discussion of the previous section, we can write the relevant field equations on 
a FLRW ansatz. We found, after some trivial algebra, that the system can be written 
as follows 

1 k 2 
a(3X - l)H + a(3A - l)H 2 = — (1 + 3w)p + Aa(l - 3A) — + Aa(3A - 1)A 2 , (12a) 

6 a 4 

k k 2 
3a(3\ - 1)H 2 + 6Aa(3A - 1)A— = p + 3Aa(3A - 1)— + 3Aa(3A - 1)A 2 , (126) 

a 2 a 4 

p + 3(w+l)Hp = , (12c) 

where we have chosen to leave the parameters a, A, A explicit because of their 
cosmological relevance, and we have defined the variable A = ^Z^x) 2 which is always 
positive, since ( < owing to the detailed balance constraint. Notice that, already at 
this level, only a subset of all the initially defined couplings (a, A, A, A), plays a role in 
the cosmology. 

3.1. Dynamical analysis 

In order to analyze the phase space of this cosmological model, let us define the variables 

n = t^jh > z = -jrr > K = 2A ^m ' c = -im ' 13 

3aH z ti A a z H 2 arri 1 

with the cosmic time M = ln[a(i)]. 

The cosmological equations (112aj) - (112cp are then equivalent to the system 

1 I ^ 

n' = W M 2 + Q(2C - 2z - 3w - 1) , (14a) 
3A — 1 
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z' = -2z 


2 + z 




(146) 


K' = K 


2C- 


2z+ x n 

3A-1 




(14c) 


a = 2c 2 -c 


\2 + 2z \ +3u 'fil 
3A - 1 


(Ud) 


with the Gauss constraint 


Q 






1 + K- 


z-C 




(15) 



where the prime denotes the derivative with respect to M . The (1151) allows to eliminate 
one of the equations. If we choose to eliminate the one for Q we obtain 



z' = z[3 + C + K - 3z - 3w(C - K + z - 1)] , (16a) 

C' = C [C + K - 32 - 3w{C - K + z - 1) - 1] , (166) 

= X [1 + C + K - 3z - 3w(C -K + z-1)]. (16c) 

This system presents three invariant submanifolds z = 0, K = 0, C = which, by 
definition, cannot be crossed by any orbit. This implies that no global attractor can 
exist in this type of HL cosmology. Also, the structure of (fl5|) reveals that the system 
is non-compact and asymptotic analysis will be required in order to complete the study 
of the phase space. Due to its dimensionality it is not easy to draw a plot of the phase 
space. Naturally, since invariant submanifolds are present, one could think to give a 
plot of orbits belonging to them, however the physical meaning of such manifold is not 
clear. For example, the invariant submanifold C = might be associated at first glance 
to flat cosmology, however from the definitions ( Ti~3l . one can see that this would be the 
case only for K = 0, so that the rest of C = has no real physical meaning. Similar 
arguments can be given for z = and K — submanifolds. As consequence only in 
orbits in the phase space bulk correspond to physical evolutions for the system and we 
have describe them without the aid of graphics 

3.2. Finite analysis. 

Let us start with the finite analysis. The finite fixed points can be found setting z', C, K' 
in (1 1 6 - ( TT6cl) to be zero and solving the corresponding algebraic equations. The results 
are shown in Table [U The solutions associated to the fixed point can be derived from 
the Raychaudhuri equation 

H = ~[3z-C -K -3 + 3w(C - K + z - 1)] H 2 , (17) 

and the results are shown in Tab. [H too. As one can see there, we have a general 
Friedmann solution which depends on the barotropic factor w of standard matter, a 

|| Of course one could consider the O = invariant submanifold which could be of interest, but in our 
setting this would require the derivation of a new system of the type (|16a| - (|16c| in which another 
variable is eliminated. Since here we are mainly interested in the role of the HL corrections, drawing 
= will not add much to the understanding of the dynamics of the cosmology. 
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pure radiation like solution, a Milne universe, and an exponential solution. Note that 
the Friedmann solution is generated when the term related to the spatial curvature k is 
dominant. Instead, when the effective Hofava radiation is dominant, the corresponding 
cosmology presents a "radiation-dominated like" solution, as expected from the form of 
the corresponding terms in (112aj) - (112cj) . Substitution into Eqs. fll2ap -( fl2cj) reveals that 
this is actually a vacuum solution and that the exponential solution corresponds to a 
an oscillating evolution with period T = 2n^j^- ^J^. Such solution can be connected 

to the scenario proposed in [12]. The Milne solution, instead, is not an actual solution 
of the syste nfl 

The stability of the fixed points can be determined by evaluating the eigenvalues 
of the Jacobian matrix associated with the system f ll6flP -( [T6cj) . as prescribed by the 
Hartman-Grobman Theorem [27] . The results can be found in Table [TJ We can observe 
that the thermodynamical properties of matter influence the stability of the Friedmann 
point A and of the point £>. This means that, if < w < 1/3, the typical completely 
finite orbit implies an initial radiation-like behavior that evolves towards a Friedmann 
or a Milne behavior (or both), to eventually approach an oscillating state. Instead, 
if 1/3 < w < 1, A is a source so that the typical orbit will start with a Friedmann 
evolution and evolve towards a radiation-like or a Milne evolution before converging to 
an oscillatory behavior. In both these scenarios we do not find any transition to a dark 
energy era. 

In fact, in terms of the phase space one can characterize this transition as the fact 
that the orbits will cross the plane 



but none of the finite orbits described above can actually achieve that. This conclusion 
was in some sense expected, since it was already noticed [11] that in this case the 
cosmological term seems to have the wrong sign when the detailed balance condition is 
imposed. 

3.3. Asymptotic analysis 

Owing to the fact that the dynamical system (13.11) is non-compact, there could be 
features in the asymptotic regime which are non trivial for the global dynamics. Thus, 
in order to complete the analysis of the phase space we will now extend our study using 
the Poincare projection method |28j . This method is based on a phase space coordinate 
change 





C = psinO cos0, 



K = psin#sin0, 



z = p cos 6 , 



(19) 



where 



r 



(20) 



P 



1 + r 



% However this does not constitute a real problem because, as we will see, this point is always unstable 
and there is no orbit which can reach it. 
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Table 1. Finite fixed points of the system (116a|) - (|16c|) and their associated solutions. 



Point Coordinates Solution Energy density Stability 

[Q, z, C, K] 



A [3A- 1,0,0,0] a = a (t-t ) 3(1+ » ) P = Po(t - t ) 



B [0,0,1,0] a = a (t-t )2 p = 



-2 



T = I 



ICI A 
a 3A-1 



saddle < w < 1/3 

repeller 1/3 < w < 1 

repeller < w < 1/3 

saddle 1/3 < w < 1 



C [0,0,0,-1] a = a (t-t ) p = saddle 

V [0, 1, 0, 0] a = a e r( '~' o) p = attractor 



9 E [0,vr], and <p e [0,2vr]. 

In this way, the limit p — > 1 corresponds to r — > oo, and one can study the 
asymptotics of the system. Performing the transformation (T191) . the system (11 6 aft 
becomes, in the limit p — > 1, 

p — > 3(1 + w) cos6* — sin^ [(1 — 3w) cos0 + (1 + 3w) sin0] , (21a) 

0' -> - sin(2^) [3 + cos(20)] , (216) 

<p' -> sin(20) . (21c) 

Notice that the radial equation does not contain the radial coordinate, so that the fixed 
points can be obtained using the angular equations only. Setting <// = and 6' = 0, we 
obtain the fixed points which are listed in Table [2j The solutions associated with these 
points can be found with the method described in [23] . The corresponding results are 
also given in Table [2j 

The stability of these points is studied by analyzing first the stability of the angular 
coordinates and then deducing, from the sign of Eq. (121 ap . the stability on the radial 
direction [23] . Ensuing results are given in Table [3j The presence of these additional 
points greatly enriches the dynamics of the cosmology. The new points represent a 
Gauss function evolution that is characterized by a very fast expansion and a similarly 
fast recollapse. If unstable, such points could represent — at least in principle — an 
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Point Coordinates Solution 

(M 

Ex [0, 0] 

S 2 [tt, it] 

Tx [0,tt/2] 

T 2 [vr, tt/2] 

Qi [tt/2, tt/2] 

£ 2 [3tt/2, tt/2] 

Table 2. Asymptotic fixed points of the system (|16a|) - (|16cp . their associated solutions, 
and stability. 

inflationary or dark energy era. However, none of the orbits connecting these points 
can be considered of much cosmological relevance, because the ones that involve the 
asymptotic fixed points with the right stability behavior do not contain a transient 
Friedmann era. 



a = a e§ («+i)(*-*o) a 
a = a ef ( w+1 «*-*°) 2 
a = aoe^ 1 - 3 ^*-*") 2 
a = ao eKi-^)('-*o) 2 

a = a ei( 1+3u ')( t - t( >) 2 
a = aoe i(i+3-)(*-*o) 2 



4. Hofava-Lifshitz cosmology: no detailed balance case (SVW potential) 

If we do not impose detailed balance to hold, the cosmological equations in presence 
of matter (assumed to be a barotropic fluid) can be written as (we closely follow the 
notation of [9]) 

1 - 3 4) (h + H*) + l^p(l + 3w) - £ + ^1 + = o , (22a) 



where 



2 J \ J 2 ry '6 6a 4 3a 6 

l-^W-^-« 2 P-^-^r-4 = 0, (226) 

2 7 6a 2 P 6 6a 4 6a 6 V ; 

p + 3(> + l)iZ> = , (22c) 



= -L , (23a) 

Xi = , 236 
o 

X2 = -Qgiu 2 > , (23c) 

X3 = 12a {3g 2 + 9s) , (23d) 
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Table 3. Stability of the asymptotic fixed points of the system (|16a|) - (|16c[) . 



Point 


Eigenvalues 




p' 


Stability 


Si 


[2iV2, -2iy/2\ 


p'> 


always 


attractive center 




\2i\/2 —2i\/2\ 


0' < 


\j c\ii vv ex y 


rpnnl^i vp ppntpr 


Ti 


[2V2, -2y/2\ < 


p > u 
[p'<0 


fl < in <r 1 
U _^ UJ \ 2 , 

\<w<\ 


saddle 


Ti 


[2V2, -2V2\ < 


\ ft >o 
[p'<0 


\ < w < 1, 
< w < §, 


saddle 


Gi 


[2i, -2i] 


p'< 


always 


repulsive center 


G-2 


[2i, -2i] 


pl> 


always 


attractive center 



X4 = 24 (904 + 3g 5 + g 6 ) , (23 e) 

and we have denned £ = 1 — A assuming also that £ 7^ 2/3. Notice that taking £ > 2/3 
would imply that the energy density in the Friedmann equation has a negative sign. In 
fact, this range of the parameter corresponds to spin zero modes of the theory, which can 
be excluded at the cosmological level and are related to unwanted ghost modes. Also 
the sign of the term \i determines the sign of an (effective) cosmological constant in 
the model and \2 can be always taken to be negative p]. Comparing the system above 
with the one in ( ll2fflP -( fT2cj) one can note that, as pointed out in [9], there is not much 
difference between the two cases. For example, apart from the values of the constants, 
(I12fej) and (122 fej) differ only by the term associated with % 4 . However, as we are going to 



show, the associated cosmological dynamics will be non-trivially changed. 

4-1. Dynamical analysis 
Let us define the variables 

«V x = Jfx^_ V = J^_ Z= *L_ K = ^- (24) 
# 2 ' 6a 4 # 2 ' y 6a 6 # 2 ' 6# 2 ' 6a 2 # 2 ' K ] 

and the cosmic time Af = \n[a(t)]. The cosmological equations (12 2 aft - (12 2 eft are then 

equivalent to the system 

= U(z- X -2y) _ _ 1 _ 2(1 + 3 W ) 

2-3£ J 2-3£ y J 
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2x 4x 2 

x ' = J^~ Ay + 2z ~( 1 + 3w )^ ~ 2 ] " ^£ . ( 25& ) 

2y r -2x + 2z - (1 + 3w)n - 4] - , (25c) 



2-3^ 2-3£ 
2z 4z 2 
z ' = 2^3£ [ ~ 2x " % " (1 + 3w)n + 2] + 2~^3£ ' (2M) 
o if 

^' = -^3^[2^ + 4y-2^+(l + 3 W )fi] , (25e) 

with the Gauss constraint 

K + x + y + z + £l+l-~Z = 0, (26) 

where the prime denotes derivation with respect to M . As before, the constraint ([TBI) 
allows to eliminate one of the equations. If we choose to eliminate the one for Q, we 
obtain 

2(1 - 3w)x 2 
2-3£ 

[2K{3w + l)-6(w-l)y-6{w + l)z-{l-3w)(2-3£)) , (27a) 



x 



2-3£ 



y 



6(w — l)y 2 



z 



2-3£ 

i 

6(w + l)z 2 



K' 



2-3£ 

— \2K 

3C 

2K 2 (3w + 1) 



+ — ^— [2K(3w + 1) + 2(1 - 3w)x - 6(w - l)y + 3(w + 1)(2 - 3f)] , (27c) 



2-3^ 

+ —^—[2(1 - 3u>)x - 6(w - l)y - 6(w + l)z + (3w + 1)(2 - 3£)]- (27d) 

As expected, the new degrees of freedom, associated with additional terms present in 
this case, result in an additional dimension for the phase space. The system above 
possesses four invariant submanifolds, namely x = 0, y = 0, z = 0, K = 0. This implies 
that, also in this case, no global attractor can exist. As before, the structure of ( 1261) is 
such that the system is non-compact and this will require an asymptotic analysis to be 
performed. Also in this case, the dimensionality of the phase space and the nature of 
the invariant submanifolds is such that we cannot picture clearly the phase space. Thus 
we will not be able to use graphics to show our results. 



4-2. Finite analysis 

The finite fixed points are found by setting x',y,z',K' = in (j27ap - (]27d|) and solving 
the resulting system of algebraic equations. The results are shown in Table [H Note 
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that the presence of a fixed point depends on the exact sign of the constant Xu as wen 
as on the value of £. For example, given the fact that the variable Q is defined positive, 
the fixed point A, can exist only if its coordinates are non negative, and this happens 
for £ < | only. 

The solutions associated to the fixed point can be derived from the Raychaudhuri 
equation 

H 2 

H = — \2K + 2x + 6y - 6z - 6w{x + y + z - 1 - K) - 9(w + 1)£ + 6] . (28) 

4 — 64 

The corresponding results are shown in Table (j4j). The new terms in the system (122 ah - 
(I22cl) induce a new fixed point characterized by a behavior t 1 ^ 3 which corresponds to the 
domination of a new cosmic component which goes like a -6 . Substitution into Eqs. (1 1 2 aj) - 
(I12cf reveals that, in this case, only the Friedmann solution and the exponential solution 
yield identities. Specifically, the exponential solution represents a de Sitter solution if 
> 0, otherwise it is associated with oscillations. In other words, if one wants 
standard matter to interact with HL gravity in the standard way (gravity makes matter 
to attract itself), then the de Sitter solution is only present if Xi > 0, i.e., if the 
cosmological constant has the right sign, as expected. 

In the same way as in the previous case, the stability of the fixed points can be 
determined by evaluating the eigenvalues of the Jacobian matrix associated with the 
system (II 6 aft , as prescribed by the Hartman-Grobman theorem [27] . The results can be 
also found in Table HI The stabilities of the fixed points in this model are different from 
the corresponding ones in the previous case although, again, the only stable finite fixed 
point continues to be the de Sitter one. 

As before we can characterize the transition form decelerated and accelerate 
expansion as the crossing of the hyperplane 

= 2K + 2x + 6y - 6z - 6w{x + y + z - 1 - K) - 3(1 + 3w)£ + 2 . (29) 

Unfortunately, due to the higher dimensionality of the phase space, the dynamics 
of this model are not as easy to extract as the ones in the previous paragraph. However 
in the general structure of the equations there is no feature which prevents the existence 
of an orbit connecting the unstable Friedmann phase with the de Sitter attractor. This 
means that HL cosmologies without detailed balance can actually admit a transition 
between Friedmann evolution and a dark energy era. 

4-3. Asymptotic analysis 

Let us now look at the asymptotic analysis of the dynamical system (125 a[) - ( T2~5el) which, 
like (13. ip . is non-compact. We will use again the Poincare projection method [28]. Let 
us perform the phase space coordinate change 

x = p cos sin ^ sin?/', y = p sin ^ sin sin?/', z = cos # sin?/', K = pcosip, (30) 
where again 
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Point Coordinates Solution Energy density Stability 

[Q,x,y,z,K] 



A 


rifo _ n n n nl 


Li 


— n „(f — 

— a 0\ L 


■1 

/ 


p — 


Po\i to) 


sdu.u.ie 


B 


[0,|(2-3O,0,0,0] 


a 


= ao(t - 




p = 





saddle 


C 


[0,0,0,0,1(3^-2)] 


a 


= a (t - 


to) 


p = 





saddle 


V 


[0,0,0,|(2-30,0] 


a 


= a e T ^- 


-to) 


p = 





attractor 


s 


[0,0,i(2-3O,0,0] 


a 


= a (t - 


to)* 


p = 





repeller 



T 



3(2-30 



Table 4. Finite fixed points of the system (|27a[) - (|27dl) . their associated solutions and 
stability. 



9 G [0, 7r], and G [0, 27r]. Performing the transformation ( l30l . the system f l27a!) -( l27d[) 
becomes now, in the limit p — > 1, 

2 

p' — > -{sin #[(1 — 3w) cos0 + 3(1 — w) sin0] — 3(w + 1) cos#} sin-0 

2 — 34 

2(3w + l) 

^cos^, (0) 

6'^ -sin(20)[cos(20) -5] , (1) 

0' - sin(20) , (2) 

if/ -> -i sin(2?/>) [2 cos 2 + (cos(20) - 3) sin 2 6] . (33) 

Also here the radial equation does not contain the radial coordinate, so that the fixed 
points can be obtained using just the angular equations. Setting 0' = 0, 0' = and 
ip' = 0, we obtain the fixed points listed, together with their solutions, in Table [51 Their 
stability behavior is given in Table El 

5. Discussion and conclusions 



In this paper we have used dynamical system techniques to analyze the non vacuum 
cosmology of Horava-Lifshitz gravity both in the presence and in the absence of detailed 
balance. Our analysis has allowed to both gain an understanding of the qualitative 
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Point 


Coordinates 

[<f>,94] 




Solution 


Ti 


[0,0,0] 


a 


= aoe^ (1+3tu)( *~* o)2 


T 2 


[7T,7T,0] 


a 




Gi 


[0,0,vr/2] 


a 


= a e§ (-+!)(^*o) 2 




[0,7T,7T/2] 


a 


= a el^ +1 )(*-*°) 2 


n x 


[7r/2,7r/2,7r/2] 


a 




Ho 


\ir/2 7r/2 37r/2l 


Q 


= a e5( u, + 1 )(*-*°) 2 


Ti 


[0,7r/2,7r/2] 


a 


= aoe^ 1 - 3 -)^) 2 




[vr,7r/2,7r/2] 


a 





Table 5. Asymptotic fixed points of the system (|27aj) - (|27rfl) and their associated 
solutions. 

behavior of the cosmology and to obtain interesting new constraints on the parameters 
of the model. 

In the first case the phase space exhibits four finite fixed points, three of which 
represent physical solutions of the system. Although one of these points is associated 
to an unstable classical Friedmann solution that could be certainly useful to model the 
nucleosynthesis and the structure formation periods, our analysis does not reveal any 
useful fixed point which could model an inflationary or dark energy phase. It has been 
proposed that, because of the changes in the value of the speed of light contained in 
the theory, the absence of an explicit inflationary phase might not be such a serious 
shortcoming of the theory [8], although at first glance it is difficult in this setting to 
produce a dark energy era. A more conclusive analysis of this issue, however, will require 
a complete numerical study, which is left to future work. 

An interesting result of our investigation is that one of the attractors in the theory 
corresponds to an oscillatory behavior. Such oscillations can be associated with a 
bouncing universe, which can be connected to the analysis in [12]. Other attractors are 
associated with a recollapsing solution in the form of a Gauss function. These solutions, 
which also appear in the context of GR plus a scalar field, higher-order and scalar- 
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Point 


Eigenvalues 




Stability 


Tx 


[2y/2, -2y/2, -2] 


p' < 


saddle 


T 2 


[2V2, -2V2, -2] 


p' > 


saddle 


Qi 


[2y/2, -2V2, 2} 


p' < 


saddle 


Go 


[2>/2, -2V2, 2] 


p' > 


saddle 


n x 


[-4,-2^3,2^3] 


p' < 


saddle 


n 2 


[-4,-2^3,2^3] 


p' > 


saddle 



2i [2i>/2, -2i>/2, -2] p' < < 



£< 



2 
3 



< w < I attractive center 



~ < w < 1 
< w < \ 



repulsive center 
repulsive center 



£ > § \ < w < 1 attractive center 



X 2 [2zv^, -2zV2, -2] p' > < 



C < § 0<w<| repulsive center 

£ < § 3 < w < 1 attractive center 

C > § 0<w<| attractive center 

£ > § !<il><1 repulsive center 



Table 6. Stability of the asymptotic fixed points of the system (|27ap - (|27rf|) - 



tensor gravities, could — if proven to be unstable — model some kind of superinflationary 
phase whose properties ought to be carefully investigated. However, most of the orbits 
that contain these points will either evolve towards recollapse or to the oscillatory 
attractor. In other words, when detailed balance is imposed, although one is able to 
prove that the corresponding cosmology could, given suitable initial conditions, behave 
like a Friedmann one, it turns out that it will eventually evolve towards an oscillatory 
behavior or either recollapse, leaving no space for a dark-energy era. 

Maybe the most important result in this paper is related to the HL cosmology 
without detailed balance. In this case, in fact, the additional freedom in the values of 
the parameters allows the existence of cosmic histories which contain a Friedmann era 
and evolve towards a dark energy one. This follows because the phase space contains a 
fixed point associated to the standard Friedmann solution which is unstable and another 
one which can be associated to a de Sitter solution which is an attractor. The last point 
can then model a dark energy era. However, the existence of these fixed points is only a 
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necessary condition: because of the presence of invariant manifolds and the constraints 
on the parameters only a subset of the phase space and the parameter space will realize 
this scenario. 

On the other hand the fact that the fixed points all lie on invariant submanifolds 
guarantee that such orbits can exist. Unfortunately the high dimension of the phase 
space makes it quite hard to perform any qualitative analysis. Therefore only numerical 
methods will allow the investigation of the details of these orbits. Notwithstanding 
these problems, we feel that it is safe to conclude that a cosmology compatible with the 
present observations can be obtained, in the HL framework, only if the detailed balance 
is broken. Such result makes this type of HL gravity a very promising phenomenological 
model for both the study of dark energy and quantum gravity. 
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